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(N ■ Abstract 

. The gauge-independent invariant approach to investigation of the linear scalar perturbations of 

ry^ ', inflaton and gravitational fields is developed in self-consistent way. This approach allows to compare 

\ various gauges used by other researchers and to find unambiguous selection criteria of physical 

and coordinate effects. We have shown that the so-called longitudinal gauge commonly used for 
studying the gravitational instability leads to overestimation of physical effects due to the presence 

. of nonphysical proper time perturbations. Equation of invariant dynamics (EID) is derived. The 

general long-wave solution of EID for an arbitrary potential U{4>) has been obtained. We have 
also found analytical solutions for all wave lengths at all stages of the universe evolution in the 
framework of simplest potential U(ip) = m 2 tp 2 /2. We have constructed the analytical expressions 
for the energy density perturbations spectrum A(k, t) at all possible k and t. Amplitude of the 
long-wave spectrum in the case of the transition from short waves to long ones occurs at the 
inflationary stage is almost flat, i.e. has the Harrison-Zeldovich form, for arbitrary potential U((p). 
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I. INTRODUCTION 



The standard inflationary model 0, 0, IH is the most successful theoretical model for 
explanation of the observable Universe structure. According to this model at the early 
stages the Universe was in an unstable vacuum-like state characterized by a slow linear 
drop of the Hubble parameter Hit) with time growth thus the cosmological expansion had 
a quasiexponential character 



a(t) ~ exp J H(t)dt. 



This stage of evolution of the Universe is called the inflationary epoch. For reviews on 
inflationary cosmology see Refs. jHIHEl]- 

The simplest model of the cosmological inflation for the flat universe is described by 
Lagrangian 

C = ~R{gik) + k<P, i - U(0). (1.1) 

zx A 

Einstein equations for this model are 

R k - ~$h = xTf , f t k = (p,^ k - 5 k (hp.^ - U(0)j , (1.2) 

where x = 1/Mp l is the gravitational constant, Mpi is the Plank mass. The equation of 
motion of the inflaton field is 



dim 

dip 



^;: + ^F = o. (i.3) 



The symbol """ here implies that metric and inflaton (scalar) field (p are quantum 
operators. We will divide their into the classical spatially homogeneous parts denoted g^ and 
(pit) correspondingly and the quantum fluctuation operators x(t, a?) and hik(t,x), describing 
small perturbations. 

The main properties of the inflationary universe may be seen in the model with quadratic 
potential 0] 

U(0) = \m 2 (p\ (1.4) 

This model is used in our work for describing of some new aspects of the self-consistent 
invariant dynamics of inflaton and gravitational fields. 

Traditionally there are three main problems in the framework of the model (jl.lj) . The 
first problem is concerned with the self-consistent dynamics of spatially homogeneous fields 

= 9ik = (gik) = diag(l, -a 2 (t), -a 2 (t), -a 2 (t)). 

that is described by equations (jl.2j) and (jOj) in the zeroth order of perturbation theory: 

3H 2 = x (h 2 + U((p) 

[1.5) 

dUiy) 
(p + 3H(p+ — ^ = 0. 

G(p 
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Here H = a/a. In further calculations we'll use both the cosmic time t (goo = 1) and the 
conformal time 77 (g o = ° 2 )- Recall, the dot " ' " denotes a derivative with respect to the 
cosmic time, the prime " ' " denotes one with respect to the conformal time. 

Solution of this problem is well-known and underlies the chaotic inflation scenario 0. 
We have unified the corresponding calculations. The system of background equations (|1.5|) 
is solved analytically at all times t. Consideration of dynamical properties for the simplest 
inflation model is necessary for solution of the second problem: 

The second problem is connected with the self-consistent dynamics of spatially inhomo- 
geneous quantum fluctuations of inflaton and gravitational fields 

X(t, x) = <p(t, x) - (p(t), h ik (t, x) = g ik (t, x) - (g ik ) 

both at inflationary and postinflationary stages of the Universe evolution. The theory of 
cosmological perturbations is based on expanding the Einstein equations to linear order 
about the background metric. The theory was initially developed in pioneering works by 
Lifshitz [7|. Two aspects of the problem are: (i) the stability of the inflationary process; (ii) 
the back reaction of the long-length perturbations on the expansion of the Universe in the 
past and in the present epoch 

The third problem is how the Harrison- Zeldovich spectrum Q was formed from the 
inflaton field vacuum fluctuations, the length of which is much less than the Universe size 
at the beginning of the inflation. As it is known, the solution of this problem is the basis 
for the modern theory of the large scale structure formation in the Universe. 

A large amount of papers is devoted to solving the last two problems mentioned above 
which has been investigated by various researches for more than twenty years. The con- 
ventional approach is based on the investigation of the scalar perturbations in the so-called 
diagonal or longitudinal gauge as the most comfortable. The theory is formulated in terms of 
the relativistic scalar potential t) = Sg 00 (x,t)] all other observable values are expressed 
via t). The subject of the discussion, which constantly appears in the literature, is 
if the predictions of the physical consequences of the theory are invariant. Two point of 
view clashed at this point. Some people suppose, that the long itudinal gauge is physically 
preferred because its main object, £), is invariant itsel f llOi . The opposite claim is that 



t) describes effects in the fixed frame of reference |11|,|12||; in other frames of reference 
the observed phenomena may look in another way. 

In fact, the question is if the theory in the longitudinal gauge can be used to reconstruct 
the past of the Universe on the basis of the observations that will be performed. The 
Fourier-image of $k(£) is the functional $k{i?(t)}, and the Hubble function H(t) contains 
information about the inflaton field and the inflationary process itself. The theoretical 
reconstruction of the past is possible only in the case of &k{H(t)} is an invariant functional. 
On the other hand, in the case of the information contained in ^^{Hft)} essentially depends 
on the properties of the frame of reference that is prescribed by the longitudinal gauge, 
the theory of relativistic scalar potential can not be used to interpret the results of the 
experiment. 

Our claim — that the theory in the longitudinal gauge is noninvariant — is made on the 
basis of the proposed invariant dynamics of the scalar perturbations. The main result of 
our investigation is the strict proof of the following statement. The invariant information 
about the dynamics of the scalar perturbations is selected from the equations of the linear 
gravitational instability theory by the identical mathematical transformations without making 
use of any gauge condition at any stage of the calculations. Our theory is formulated as a 
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closed system of equations for the invariant Jk of metric perturbations, invariant function 
Xkinv of the perturbations of the inflanton field, its derivative Xkinv and energy density 
perturbations 5e- mv . The theory is such that the invariance of the physical values follows 
from their mathematical definitions, and the invariance of the equations does from the way 
to obtain them without any gauge. 

Invariant approach allows us to compare different gauges which are used in the works 
of other researches and to find unambiguous separation criteria of physical and coordinate 
effects. The problem of such criteria existence was widely discussed also in j^, [1(1 [HI, ■ 
We have shown that the longitudinal gauge leads to the overestimation of the physical effect 
as a result of the strong perturbations of the proper time in frame of reference specified 
by the longitudinal gauge. The general qualitative properties — the power-law instability 
of long-length perturbations and the formation of the Harrison-Zeldovich spectrum — are 
the same in the two approaches, but the numerical values are different, the perturbations 
in the longitudinal gauge theory being several times greater. Invariant approach excludes 
nonphysical coordinate effects and gives a key for analytical investigation of equations in 
the perturbation theory that is the aim of our work. 

Among the theories formulated in the fixed gauges, the synchronous gauge theory pro- 
posed by Lifshitz 0] is an adequate one, because it this gauge the evolution of the background 
and perturbations is analyzed in one and the same time. 

We have analytically investigated of the equation of invariant dynamics (EID) for Jk. 
The general long-wave solution of EID as a functional of background solution for arbitrary 
potential U (0) is obtained. For one of the simplest model of inflaton field (|1.4|) we have found 
explicit solutions of EID for all wave numbers k and times t. An important point of our 
investigation is that both the background characteristics and characteristics of perturbations 
at the inflationary and post inflationary stages in the considering model (jl.lj) of the Universe 
are completely defined by the parameter m and the initial values of the Hubble function Hq 
and the scale factor So- 
il. BACKGROUND DYNAMICS: AN ANALYTICAL GLANCE AT THE LINDE'S 
CHAOTIC INFLATION 

The system (|1.5|) for potential (jl.4|) turns to 




(2.1) 



(p + 3H(p + m 2 ip = 0. 
In the limit t — > solutions of this system can be represented as power series 



CO 



oo 





(2.2) 



n=0 



ra=0 



Substituting ()2.2|) in the (j2.1|) we obtain the system of equations 



oo 



oo 




n=0 



k,n=0 



(2.3) 



oo 



CO 



OO 



n(n - l)a n t n - 2 + 3 ^ na n b k t k + n ~ x + m 2 



J2a n t n = 0. 




n=0 
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The initial value of H(t) unambiguously defines the initial values of functions tp(t) and 
<p(t). For convergence of series we should assume that 



2 

m 
96f 



< 1, m < 1, a 2 = 0. (2.4) 



Taking into account these conditions we will receive from (|2.Hjl the following relations for a r , 
and b n : 

u / a \ V2 / a \ V 2 5 / a \ V 2 o 



a °=mUJ ' ai = _ yUJ ' a3= 162l|UJ ' a « = --^n-i for n^4; 

&! = , 6 2 = 0, b n = —xaia n for n ^ 3. 

3 

Therefore the expression for expansion coefficient of function H(t) is 



6m6 (36q 



.71 — 3 



816§ n! 



The series can be summed up and the result is 



, , m 2 m 6 t 2 2m 6 , n h , \ , 

= <*> - - ( + sW " m ( e ' 1 + 3M > ■ (2 ' 5) 



For small £ due to conditions ()2.4|) this expansion corresponds to the linear approximation 
- peculiar feature of the inflationary epoch (Fig. CJ: 

2 

771 

H(t) = H - — t. (2.6) 

Let us build an analytical solution of (|2.1j) in the limit t — > oo, corresponding to the 
postinflationary stage. As it is known in this limit the Hubble function H(t) must approach 
to 2/3t of present-day Friedman's universe. In the equation Cp + 3Hp> + m 2 ip = one can 
make the substitution p = ip/a 3 ^ 2 . As a result we have 

^ + ^{m 2 -^-H-^H 2 ^ =0. (2.7) 

This is the equation of oscillator with the variable frequency u>(t) 2 = m 2 + p(t), where 
pit) = —\H — \ H 2 — > steadily for t — > oo. The general solution of equations like ()2.7|) 
has the form: 

The function e{t) should be defined. Substituting ()2.8|) in (|2.7j) and differentiating we obtain: 

j t {m 2 + I)D = 0, (2.9) 
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where D(t) = l/e(t), I is the integro-differential operator, which is defined as 
for example, II — p/2, Ip — \{'p + 3p 2 ). Equation ()2.9|) gives 



Constant has been chosen for correspondence to zeroth approximation of equation ()2.7|) - 
the harmonic oscillator with the frequency m. Suffice it to take into account just a first 
correction. Finally we have 



mil- -^—(6H + 9H 2 ) ) . (2.10) 



I ! i / / — i- ' J / / 

im? 

Then instead of <p we introduce the new function G = ip 2 , as a result we get 

H (\ 



3H 2 



a 3 



{^G + m 2 G-\{HG)^. (2.11) 



On the other hand using (|2.8j) and (J2.1(Jj) we have up to the first order terms of 1/t inclusively 



1 A /„ . . x (ZH , 9 



G = ip 2 = — [1 + cos(2mt + a) + — + — / H 2 dt sin(2mt + a) ) . (2.12) 

Substituting this expression in (j2.11|) we get an equation for H(t): 

k ( 3H 



3H 2 = ^- i m + — sin(2mt + a) ) . 
a 3 \ 2 J 

This equation has been solved by substitution a(i) = e' H ( l ) dt then by taking the logarithm 
and differentiation. Up to the second order terms of 1/t inclusively we get the following 
asymptotic solution of the system (j2.1j) for t — > oo: 

H(t) = + M^F in(2m(t - io) + a) - (2 - 13) 

This solution corresponds to the postinflationary stage of the Universe (Fig.^). Two param- 
eters — shift parameter to and initial phase a — can be found by matching of background 
solutions J2IED and (f2~T3j) . 

Let ti is the transition time of the Universe from inflationary to postinflationary stage. 
At t — t\ we require an execution of the following three conditions: 

1. Continuity condition of H(t): 

m 2 2 1 

H Q - — ti = -r- — + — — — g sin(2m(t! - t Q ) + a); 

3 3(ti — to) 3m(ti — to) 
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FIG. 1: Hubble function as a function of cosmic time at the inflationary (left) and postinflationary 
(right) stages; Hq = 1, m = 0.15 (system of units h = c = 8ttG = 1), t\ ~ 127, to ~ 113. 



2. Smoothness condition of H(t): 
m 2 2 



3 3(t!-t ) 2 Kti-to) 2 

3. Continuity condition for envelope of H(t): 

m 2 2 

Hft — tl = — ; 

U 3 3(ti-t 



+ —j- cos(2m(ti - t ) + «); 



The solution of this system is 

t! = —^-- t Q = ^4--, a = 7r(2fc + l)-4, keZ. (2.14) 
m z m mr m 

In the model under consideration with m = 0.15 in the Planck units h = c = 8nG = 1 the 
slow-rolling regime holds form t = till t\ = 127 (Fig. [TJ. 

For subsequent calculations we need in accurate expressions for scale factor at all stages: 

1. Inflationary stage < t ^ t±: 

a(t) = a e 2 = a e 2 , (3 = —, 

2 (2 - 15) 
(#0 + 

where ao is the initial value of scale factor a(i). We suppose that H(rj) ~ const in all 
differentiations and integrations at the inflationary stage. 

2. Postinflationary stage t > t\. 

a(t) w A(t - 1 ) 7 , a{v)~-?r'n , where A = a e^ — ) . (2.16) 

9 V 2 / 
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Obviously the non-perturbed characteristics of the inflationary and postinflationary 
stages are completely defined by parameter m and initial values of Hubble function Ho 
and scale factor a - 

In inflationary scenario density perturbations are generated from vacuum fluctuations of 
the inflaton field, for reviews see Refs. In the following section we illustrate the 



main idea and obtain the basic formulas for the invariants of the metric and inflaton field 
perturbations. 



III. COSMOLOGICAL THEORY OF SCALAR PERTURBATIONS: EQUATION 
OF INVARIANT DYNAMICS. 

Proceed to the analysis of the spatially inhomogeneous quantum fluctuations: 

x{t, x) = 0{t, x) - (p(t), h ik (t, x) = g ik {t, x) - (g ik ). 

Notice, the inflaton waves are the direct consequence of that the inflaton field is described 
by nonlinear Klein-Gordon equation. The inflaton itself is usually called the spatially ho- 
mogeneous mode of this field. The spatially nonhomogeneous modes will be called below 
the inflaton fluctuations. In the context of scalar field matter, the quantum theory of cos- 
mological fluctuations was developed by Mukhanov 3 ■ 



The physical effect under discussion is in the following. According to the Einstein equa- 
tions, the fluctuations of inflaton field inevitably produce the potential fluctuations of grav- 
itational field. As a result, in the system there appears a collective motion where the 
amplitudes and phases of inflaton and gravitational wave excitations are uniquely related. 
The coefficient in the relation is the rate of spatially homogeneous inflaton change. This 
type of motion will be called below the scalar inflaton-gravitational (SIG) waves. Our work 
is devoted to studying SIG waves in the frame of the simplest model with the Lagrangian 

dm. 

It is well-known, that from the scalar perturbations of metric |K 

Sg ik = h ik = a\ V ) ( 2^"+^) ) (3-1) 



the following invariant functions can be formed [16|: 

$ = + i[(5 -£>]', 

a , (3.2) 

* = ^- -(B-E 1 ). 
a 

Classical problem of the linear perturbation theory ultimately consists in search of these 
invariants. This problem is extremely difficult and it hasn't been solved generally to the 
present day. Many authors resort to various simplifications namely the co-ordinates may 
be chosen so that the initial expression for metric perturbations (j3-H) has got more conve- 
nient form for investigations. It is reached by using of various gauges namely longitudinal, 
synchronous, co-moving gauges. But it is necessary to understand that within any gauge 
there is some arbitrariness in selecting of the co-ordinates. Therefore it is impossible to fix 
co-ordinates hard within a gauge and as a result we have not a physical effect but effects 



8 



are concerned with the co-ordinates motion itself that is the coordinates effects. Moreover 
in the investigation within any gauge it is impossible to determine the unambiguous criteria 
of the physical and coordinate effects separation. Such criteria exist only in the framework 
of invariant approach which is being developed in our work. 

The essence of our approach involves the following. Search of two invariants (|3.2j) in the 
general case for metric perturbations (|3.1|) is the despairing problem. We have seen that the 
single invariant J can be constructed from the invariants $ and 

1 fa' V I fa 2 \' 



J = $ + - -tf = (j> + - — ip) . (3.3) 
a \a J a \a J 

We will show below that for this invariant there can be obtained and solved analytically 
(exactly or by using asymptotic methods) the single second order differential equation for 
arbitrary potential U((p) without using any gauge. So, the invariant dynamics of the SIG 
waves does exist. 

We employ the notations, initially used by Lifshitz [7]. The Fourier harmonics of the 
potential excitations of metric are 

Kiy)= I i,^ U/J Kk\ | . (3.4) 

-tk p a k -(/i k + A k )5£ p~ Ak 

Operations with the spatial indices are made in all calculations with the help of 3-Euclidean 
metric. One can easily see the one-to-one correspondence between the values from (|3.1|) and 
()3.4j) . In the longitudinal gauge, which is used widely nowadays, Ok = Ak = 0; in the Lifshitz 
or synchronous gauge Qk = o"k = 0. Further when we will construct the invariant dynamics, 
we do not use neither of the these gauges. However, the results of invariant dynamics will 
show that the synchronous gauge has some physical advances, in contrast to the longitudinal 
one. 

In the first (linear) approximation the equations describe the relation between the fluc- 
tuations of metric and the fluctuations of the inflaton field: 

= \ + ft/ - hf, - 4) - h\R\ - \s\ (/$ - hj - h} m Rf) = xST? , (3.5) 



2 mr ' ir ^ \dcp 



Notice, in (|3.5j) there is no any restriction on the wave length of the fluctuations. 
From (J3.5)) we can write all equations of linear perturbation theory: 



SR° {k) - SR{k) = x5T °(k), 



SR° {k) - SR{k) ■■ 1 



2 



-A; 2 iVk + -M k - 3— Q k 
3 a a z 



k2 -N k + HL k - 3H 2 Q k , ( 3 6) 



3a 
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5R°Jk) = x£T°(k) 



3 a 



a V 3 



^ T a( k ) = — ^Xk = ^Xk, 



(3.7) 



a 



5R"M - IsZSRQt) = xST* (k) 



1 



^(k) - -<g*J2(k) = — ^ 



6a 2 



-N£ - 2-iV k + k 2 (N k + 3Q k ) + 3M k + 

(X 



+6-(M k - g k ) - 6 2 



a 

2^2 "fc 2 " 



a" a' 2 



a ar 



+ 



< + 2-N k - ^(iV k + 3Q k ) - M k - 2-M k 
a 3 a 
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(3.8) 



-iV k - 3#iV k + — (iV k + 3Q k ) + 3L k + 9HL k - 6HQ k - 6(2H + 3H 2 )Q k 



+ - 



1 k hi 3 



2 fc 2 



iV k + 3HN k - — (iV k + 3Q k ) - L k - 3#L k 



+ 



(k) = -SP a 



^ ( -^Qk^' + xL) - ( tt 1 ) Xk 



<9tA 



:Qk^ 2 + 0Xk - 



/at/ 



Xk 



Here 



N k = + A k , M k = /i k + 2ika k , L k = — = /i k + 2^^. 

a a 



(3.9) 



Equations (|3.fij) and ()3.7j) enable to express the values in T£(k) via the metric perturbations. 
In terms of conformal time the value (dU/ d(p)o can be expressed from the equation of motion 
for inflaton: 

1 

77 



1 (?£. 

<p' \d<pj 



in" a' 

^7 + 2- 

ip a 



After employing this relation and constraints (|3.fij) . ()3.7j) . equation ()3.8j) reads 



6 a 



-N k ' - 2-iV k + k 2 (N k + 3Q k ) + 3M k + 



+6-(M k - Q k ) - 6 2 



a" a' 2 



a a 



Qk+ 



4^ + ^-3^-2^ + 2^ 
3 a cr \ <p' a J \ 



(3.10) 



+ 



+ tc 



1 



2 k 2 



K + 2-iV k - ^(JV k ) + 3Q k ) - M k - 2-M k 



0. 
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The projection of these equations on the tensor basis gives two equations for the scalar 
functions iV k , M k and Q k . 

K' + 2-N k - ^(iV k + 3Q k ) - M k - 2-M k = 0, (3.11) 

do (1 



K 



- + ^)K + k 2 N k -3-Q' k -6(^--- + ^)Q k + 3-M k = 0.(3.12) 
ip' J a \ip' a a "- 1 



cp" a' 



,'2 



.a 
a 



Equation (|3.12j) is written in the form demonstrating the idea of the following calculations. 
From Eq. ()3.12|) M k may be expressed via iV k and Q k and substituted in (j3.11|) . It is the 
way to obtain the single equation for iV k and Q k . 



N{' 



k r 2-iV k - ^(iV k + 3Q k ) - M k {iV k , Q k } - 2-M k {iV k , Q k } = 0. 



a 



a 



It has a specific property: two functions iV k and Qk in this equation can be combined into 
one invariant combination only: 



Jv 



-3Q P 



(3.13) 



Coefficients of all non-invariant terms vanish via the background equations for a(rj), (p(ij). 
As a result we obtain the single second order differential equation for the invariant J k - the 
equation of invariant dynamics (EID), which in terms of the conformal time reads 



+ 



,/// 



,/2 



a" tp" 



kr + 3— - 3 2— + 2— - 2^— + 2^— - 2— 



a' 



,12 



/2 



a! p' 



(3.14) 



A = 



Let's obtain the EID in terms of cosmic time t. The value dll/d(p is suitable to express 
via the Hubble function by using background system 1)1.50 : 



1 /<9?7 



1 [/ 



1/7 



- — = = r- - 3H. 

<p \d(pj (pep 2H 

Using this result together with the equations l)3.fij) and (J3.7j) . one can obtain from 1)3.8)) the 
following equation: 



6 a 



1.2 . . 

iV k - 3#iV k + — (iV k + 3Q k ) + 3L k + 9HL k - HQ k - Q(2H + 3H 2 )Q k + 



lk 



H 



H 



+ 77-riV k + HL k + \3H Z + Q k - \2H + WV, 



3 a 



H 



3H 



+ 



(3.15) 



+ 



2~P~ 



iV k + 3HN k --(N k + 3Q k ) - L k - 3HL k 



0. 
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The projection on tensor basis gives 



k 2 

iV k + 3HN k - — (iV k + 3Q k ) - L k - 3HL k = 0, (3.16) 



a 2 



N k - ( 3H + 5r j N k + ^-N k - 3HQ k - 3 \ 2H - ™ j Q k + 3HL k = 0. (3.17) 

Of course, equations ()3.16j) . ()3.17j) turn to ()3.11|) . ()3.12|) by time transformation and use of 
background equations. 

Now let us introduce the invariant in the following form 

J k=(^)-3Q k . (3.18) 



So we find 

iV k = HJ k + ^iV k + 3HQ k , 

iV k = H j k + 2#J k + 6HQ k + 3HQ k + — JV k . 

H 

Substituting ()3.19j) into ()3.16j) and ()3.17j) . one obtains 

HJ k + {2H + 3H 2 )J k + {6H + 9H 2 )Q k + 3HQ k + ( 3H + ) N k - 



(3.19) 



H 



\%{N k + 3Q k ) - L k - 3HL k = 0, 
3 cr 



(3.20) 



3L k = J k + J k ( 2^ - 3H - § ] - 3^iV k - 9#Q k + (3.21) 
H H J H a z H 

Differentiating ()3.21|) and substituting L k { J k , iV k , Q k }, L k { J k , N k , Q k } in (|3.20|) . we will get 
the equation of invariant dynamics in the cosmic time: 

J k + | 3# + 2^-| I J k + 

(3 22) 

/V • H H 2 H H 2 '"" 
+ — + 6# + 2- - 2— - ±L + ^- - 3^- J k = 
ya 2 H H 2 H H 2 H 

This equation may also be obtained from (|3.14j) by the time transformation and using 
background equations for ip and a. 

Equations ()3.22|) and ()3.14|) are valid for any potential U(ip); a fixed potential makes 
function H(t) and its derivatives also fixed. The relations which connect the invariant 
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characteristics of the gravitational and inflaton field perturbations are 



Xkin 

H 

30 



Sep 



■h 



H 1H 



H 
2H 



J^dt 



(3.23) 



In particular, these equations are used to specify the initial conditions for invariant and its 
derivative via those for the perturbations of inflaton field and its derivative. 



IV. USE OF THE GAUGES: LONGITUDINAL VS. SYNCHRONOUS 



Now let us discuss the obtained results. One can naively to assume that the existence of 
the EID does not restrain us from using any gauge, including the longitudinal one. This, 
however, is not correct. Let us consider the relative perturbation of energy density. From 
(13. 6|) and ()3.18|) one can easily obtain 



Se 8 e^. 



+ 



5e r 



where 



(4.1) 



9H 




to 



H 



Qkdt 



-St, 



to 



where e = 3H 2 / x is the background energy density, 5r = 5 J ^fgoodt = const + \ J Qkdt is 
the perturbation of the proper time. Analogously the noninvariant parts of the perturbations 
of inflaton field and its derivative are 



»k(t), 



Xk, 



#r k (0 + \ 



Evidently the term Se non i nv does not have any physical sense. This term reflects that if 
hoo ^ then the time is not synchronized in different parts of the Universe and therefore 
measurements of the background energy density in different points of the Universe lead us to 
the different results that is the term Se non i nv describes not the physical change of the energy 
density but the change of the clock run or time flow rate. To escape of such nonphysical 
perturbations it is necessary to analyze the background dynamics and perturbation dynamics 
in one and the same time, in one and the same clock. So, it is necessary to synchronize 
the clock in the whole Universe, i.e. to put Qk = 0, i.e. to use the Lifshitz's (synchronous) 
gauge. 

Now let us analyze physical validity of using a gauge in the framework of invariant 
approach. As example we will examine the widely used longitudinal gauge. From metric 
perturbations in the general form ()3.1|) one can proceed to the metric perturbations in 
the longitudinal gauge performing the special coordinate transformations fj — r\ — B + E', 
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x % = x l + E' % . Then in the new co-ordinates the scalar metric perturbations = <\> + — (S — 
£?') + (-B — £?')', if) = ^ — — (S — £") are invariants due to (|3.2|) . Metric perturbation tensor 
in that case is diagonal: 

hm=( 2 f °iJ (4.2) 



-2if>6% 

Taking into account (|3.9|) . the equation ()3.16|) leads to equality: if) = <j> = $ k , that is we 
have a theory containing the single invariant which is named the relativistic potential. It is 
presented as an advantage of the longitudinal gauge. However, this is an illusory advantage. 
In the framework of invariant approach we have established its baselessness. 

Thus, cr k = A k = in the longitudinal gauge, at that Q k = 2$ k and /i k = — 6$ k . The 
equation for $ k may be obtained from (|3.17jl : 

' k 2 ■ HH 
a 2 H 



$k + $k I H - - ) + $ k ( - + 2H - — ) = 0. (4.3) 



Invariant (|3.18|) in the longitudinal gauge is 

Jk = ~ 6 (if) " 6$k ' (4 ' 4) 

The expression for full relative energy density perturbation (j4.1J) taking into account ()4.3|) 
and ()4.4|) has a form 

^E-inv ^E-noninv 2 k 2$ k 

7 = — + ~~r~ = ~sw $k ~ ir ~ 2$k ' (4 - 5) 

where the non-invariant part of perturbations is 

t 

^= = 2-j $ k dt. (4.6) 
e H J 

to 

One can show that the existence of this term in (J4.5|) is connected with an arbitrariness in 
the co-ordinates selection within the gauge. 

At first we consider the long- wave limit. Suppose in the longitudinal gauge we the co- 
ordinates with metric ds 2 = (1 + 2@ k )c?t 2 — a 2 (t)(l — 2Q]^5ijdx l dx^ . Here the metric pertur- 
bations are described by scalar function G k . The expression for energy density perturbations 
is 

Now let us make some co-ordinates transformation preserving the gauge [ill ]: 

t = t-\(t), x i = (l + f3)x\ (3 = const. (4.8) 

Then 

ds 2 = (1 + 26 k + 2X)dt 2 - a 2 (t)(l - 20 k -2(3 + 2HX)S ij dx i dx j , 
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so cf) = @k + A and ip — Ok + (3 — //A are metric perturbations in the new co-ordinates. 
Requirement of 0k =ip — <p gives the following expression for function X(t): 

X(t) = - ( 13 I adt + 7 J , 7 = const. 



Therefore 



a 



e k = e k + 5e k = e k - — ( > I a<it + ~ ) >. 



Then instead of ©k we introduce a new function Gk into the right hand side of expression 
(|4.7|) . As a result we have 

Obviously the expression for energy density perturbations in the new co-ordinates differs 
from corresponding one in the old co-ordinates by the value 

t 

A = 2^ J 5Q k dt. 

to 

Therefore we see that the expression 

Se inv 20k „„ „H 



t 



26k -2- / e k 



e E H 

to 

is an invariant for any coordinate transformations of the form 

Our calculations clearly prove that the use of an arbitrariness in the co-ordinates selection 
for the interpretation of physical effects is incorrect. In the longitudinal gauge together 
with real physical effects there are the effects, caused by motion of the frame of reference 
within the gauge. Selection and removing of such nonphysical coordinate effects is possible 
only in the framework of invariant approach. At Fig. |21 one has shown the qualitative 
comparison of long-wave energy density perturbations calculated at the inflationary stage in 
the both approaches. We have used the simplest model with the potential ()1.4|) in the case 
m = 0.15, Hq = 1, do = 1. Analytical expressions for 8ei nv /e are shown in Appendix. 

But if in principle we can analyze the long-wave metric perturbations in the longitudinal 
gauge, because we know a form of nonphysical terms ([4.6)1 . then the similar analysis of the 
short-wave perturbations in this gauge is absolutely impossible. Let's consider this problem 
in detail. 

Invariant theory gives the single expression for nonphysical terms (|4.6|) . which is correct 
for all wave lengths. Therefore in the longitudinal gauge the expression 



3H 2 a 2 H H 



to 
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FIG. 2: Long-wave perturbations of energy density calculated at the inflationary stage in the both 
approaches - in the invariant theory (I) and by using the longitudinal gauge (2). 

is an invariant thus it should be correct the following expression 

2 ^ k _^_ M$k _ 2 " L^o, (4.11) 



3H 2 a 2 H H 

to 

where <5$k is the metric change in transition to the new frame of reference with the gauge 
conservation. At the same time a function <5<3>k has to satisfy the equation (J4.3|) . We have 
shown that the last condition does not hold for short waves. 

Let us denote <5$k = Q and consider the relation (J4.11j) as equation for Q 

t 

Ik 2 f 

Q + -Q + HQ + H / Qdt = 0. 

3H a 2 J 

to 

We divide it now by H, differentiate, and then multiply by H. As a result we have the 
equation 

fl + tils - i) + a - **) + h± = 0, 

\ H ) \ H ) dt \3HH a 2 J 

which differs from equation for metric perturbations in the longitudinal gauge (J4.3)) by the 
term containing k 2 /a 2 . What is the reason? Probably the longitudinal gauge is restricted 
itself and an important information inevitably is lost in transition to this gauge. 

Results of SIG-waves investigations fulfilled in the various gauges are substantially dif- 
ferent. For example, in the synchronous gauge we have [l2^ : 

ds 2 = —a 2 [drj 2 — (5 a p + h a fi)dx a dx }, 

(4.12) 

#00 = ~a 2 , 9oa = 0, g a p = a 2 [(l + hG)S a0 + hik 2 G >a ^}, 
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where G is the scalar function defined in the 3-space. It satisfies equation + k 2 G = 
0. Since h° = g k = 0, the non-invariant part of full perturbations of energy density 
(14. 1J) is equal to zero for any co-ordinates transformations preserving the gauge and thus 
the calculation of 5e/e gives a proper result in the synchronous gauge in contrast to the 
longitudinal one. 

The closed system of equations for metric perturbations h and hi is obtained by Grischuk 
12l | . The equation for h is the third order differential equation and it contains the nonphys- 
ical solution h ~ H, which may be obtained by the coordinate transformations preserving 
the gauge (j4.12|) . Note that this solution is trivial for the same equation but written in terms 
of invariant ([3.18)1 . and it is easily to show that this equation coincides with the equation of 
invariant dynamics (|3.22j) . 

Now let us to say some words about the quantum perturbation theory as it is. First of all, 
notice, there is a formal problem in the theory: what function is the object of quantization. 
The pose of this problem is induced by our understanding that the coordinate effects that 
can be eliminated by the appropriate choice of a classical frame of reference must not to be 
subject to quantization. In the longitudinal gauge theory the discussion of this problem is 
again reduced to the question whether the description in terms of the relativistic potential 
is invariant or not. One of the achievements of our invariant theory is that this question is 
uniquely solved. Only the invariant metric function which is uniquely connected with the 
invariant characteristics of the inflaton field is the subject of quantization. For the short- 
wave-length metric fluctuations the normalization of quantum operator is easy to find: 



30 
~H 



2ka 



Ck exp 



k dt 



+ c_ k exp 



k dt 



(4.13) 



where Ck and c^ k are annihilation and creation operators in quantum field theory. Notice, 
expressions ()4.13|) are used only to specify the initial conditions; the quantum dynamics 
itself is described by the exact operator equation ([3.22)1 . The subject of calculations is the 
value of energy density fluctuations averaged over the Heisenberg vacuum specified in the 
beginning of the inflation: 



be 



k inv 




(4.14) 



In order to calculate the value ([4 .14)1 one can solve equation (13.22)1 under |ck| = |c^ k | = 1/2. 
The averaging over the phases of the complex numbers Ck, c_k and taking the absolute value 
corresponds to the procedure of quantum averaging (see Appendix). 



V. ANALYTICAL INVESTIGATION OF THE EID. 

At the postinflationary stage the Hubble function H(t) oscillates and there are points 
where H = 0. Thus formally the EID ([3.22)1 has to have the mathematical singularities that 
could be an obstacle in investigations both analytical and numerical. However the trans- 
formations of EID coefficients together with background equations for arbitrary potential 
p. 5)1 have shown an absence of the mathematical singularities in reality that has opened the 
possibility for analytical investigations of the EID. This internal property of the invariant 
dynamics is its significant advantage. 
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Let us introduce a new function ^k(t), related with invariant Jk(i) by relation: 



A=^* k . (5.1J 



Equation for ^k(^) has a form: 



. k + W ™ k+ |^ + ^-if + 2f-2^-§™ + 3 fl |. k = 0. (5.2) 

This equation is more convenient for analysis than initial ()H.22j) . since all terms with singu- 
larities are concentrated in the single coefficient by \I/k- 

Now let us consider the elimination mechanism of these singularities in detail. The 
background system (|1.5|) for arbitrary potential can be written in the form 

dip 

First equation has to be differentiated with respect to t twice, the second one - once, then 
obtained expression has to be multiplied by xfi. As a result we have the following necessary 
relations: 

H = --^V 2 , H = —*c(p<p, H= -xipCp 2 - xip <p\ 

KLp <p +3x0 2 ff + 3xy?y5ff + xyg 2 <9 = 0. 

dip 2 

Expressing from first three relations the all derivatives of function ip through derivatives 
of H and substituting their to the forth one, after division by 2ff we obtain the following 
identity: 

IH 2 Iff ^ 3ffff_ g 2 C%) 

4ff2 2ff 6H 2 H ~ W ' { ' 

Evidently, the structure of background equations allows to combine the terms with singu- 
larities into the smooth function that is the singularities cancel. Taking into account ()5.3j) 
the EID gains the simple form: 

*x + + (5 + ^ + 2 f - ^ + m I * k = °" (5 - 4) 

One can rewrite its in the form: 

^ k + 3ff^ k + ( % + U" + p(t)] ^ k = 0, (5.5) 



where = 0, p(t) = | (2| + 6ff) . 
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Let us introduce the characteristic function 



g(M) a\U'^{t)) + p{t)Y 

Function q(k, t) has a substantially different form at various stages of evolution but at the 
all applicable domain t 6 (0 . . oo) it decreases with the time growth. There are three types 
of metric perturbations depending on value q(k,t): 

1. Asymptotically long- wave perturbations, if 

g(M)<l; (5.6) 

2. Perturbations of middle wave lengths, if: 

g(M)~l; (5-7) 

3. Asymptotically short-wave perturbations, if 

g(M)>l- (5.8) 

Consider the long- wave limit ()5.6|) . The equation of invariant dynamics ()5.4|) in this limit 
reads 

* k + 3#* k + f ^ + 2^-2^ + 6^ ) * k = 0. (5.9) 



Taking into account (|5.3|) it is easily to show that ^k(t) = y —H / H is the particular solution 
of this equation. Therefore the general solution of the EID for long waves with the arbitrary 
potential U ((f) has a form: 

Mt) ^ + C ^Jl^ d , ,5,0, 

Deeper analysis is possible only in the framework of fixed potential. Further calculations 
are fulfilled for the simplest inflation model with the potential (|1.4J1 . In this model func- 
tion \l/k for long-wave metric perturbations at the postinflationary stage oscillates with the 
stationary amplitude: 

* k (t) ~ sin(m(i-t ) + ar/2), (5.11) 

where a is defined in (I2.14j) . This result may be also obtained by direct solution of the 
EID which in that case reduces to the Matieu-like equation where there is a time-dependent 
coefficient by the oscillating correction in the frequency. Within that context the form 
of relation (|5.11|) means the presence of the parametric resonance for long-wave metric 
perturbations at the postinflationary stage. 

Let's return to initial equation ()5.5|) and proceed to the conformal time rj: 

tt£ + 2-^ + (k 2 + a 2 U^ + a 2 p(^))*k = 0. (5.12) 
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By substitution \l/ k = £\I/ k , where £ is a new unknown function, \I/ k = V —H / H is the 
particular long wave solution, this equation is reduced to 

e + 2e ^f + ^ = °" (5 - 13) 
At the inflationary stage this equation can be reduced to the Bessel equation for any k: 

_ 2 v — ^__L^ + ^ = 

77 



Its exact solution is 



m 2 1 , \ / m 2 I 



Ci BesselY ( 1 + — + -, krj J + C 2 BesselJ I 1 + — + -, krj 



At the most part of inflationary stage the inequality <C 1 is valid. In this case we have 

£(77) = Li(kr) cos(krj) — sm(kr))) + L 2 (kri sin(kr)) + cos(kr))). 
Therefore the general solution of the EID at the inflationary stage has a form 



\T-H ( f k , / fk ,\ . f f k 



».(*)=£.— / "* cos -sin /-«]] + 



r { f k . . ( fk , \ f f k . 

+L 2 — 77 — I / —at sin ( / —at | + cos | / —at 



(5.14) 



if \.l a \J a J \J a 

Entry conditions for \I / k(^) 




x 1 • xk 1 



^ V 2A;5 ' * k0 V 2 5 2 
give the following expressions for integration constants L\ and L2: 

2 /3>{ cos(krjo) — sm(krjo) 2 /3xcos(/c77 ) + sin(A;77o) 




;^ ■ ^=-^vt ;^ ■ (5 - i5) 



where 



% = T](t = 0) = J 



dt 
a 



1 



r^~> 

t=0 CLqHq 



The equation q(k, t) — 1 represents relation between /c and t at the boundary that sepa- 
rates the evolutional stages. Its solution t(k) is the time of transition from short waves to 
long ones for every k. There are two variants depending on k: 

1. The transition from short-wave perturbations to long- wave ones occurs at the infla- 
tionary stage that is < t(k) < ti; 

2. The transition from short-wave perturbations to long- wave ones occurs at the postin- 
flationary stage that is t(k) > ti- 
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For every variant the EID (|5.4|) at the postinflationary stage should be solved separately. 

In the first case at the times t > t\ there are long-wave metric perturbations described by 
()5.10|) or ()5.11|) . After matching with solution ()5.14|) at the moment t\ when the inflation 
terminates, we obtain 



Mt) = -Hi 



3>c cos(kr] ) + sin(A;?7o) V~H 




mk 3 / 2 



H 



(5.16) 



Finally consider the second particular case. The short-wave solution at the inflationary 
stage is ()5.14j) . Let us obtain the short-wave solution at the postinflationary stage. Let us 
introduce a new function (ft = a\l/k into the ()5.12j) . Using ()2.16j) . we get 



+ 



, n 2 m 2 A 6 

At H 

2 81 



4mA 3 



-1] sin 



/ 2mA 3 



27 



rf + f3 



rj* 81 3 
In the short-wave limit on right t — > t\ this equation can be simplified: 

,2 



0. 



(5.17) 



+ 6 [h* + 



m 



0. 



(5.18) 



Obviously this equation is an oscillator-like equation and it should be solved by the asymp- 
totic method which has described in the Section 2. One has obtained the following general 
asymptotic solution: 



*k(*) 



1 



(t ~ to) 2 / 3 yfl + *$-{t - t ) 4 / 3 - - to) 8 / 3 



Mi cos ^3pm 1/3 (t - t ) 1/3 + 



3m 5/3 

lbp 



(t-t ) 5/3 -^(t-t 



m 



(5.19) 



M, sin ( Spm^t - to) 1 / 3 + ^(t - t ) 5/3 - £{t - to) 3 

lOp 8p d 



where p = p(fc), 



Mi 



M 2 



2x 




18^o 2 



m 2 / 3 V P ag/ 2 m 5 / 2 exp 



2x 




cos 



COS 




+ sin 



sm 




(5.20) 



The equation (|5.4jl for function x = \I / k(t — to) has a form: 

x "+(l + ^-isin(2T+a))x = 0. 



(5.21) 



Here t = m(t — 1 ) is the new dimensionless variable. Depending on r part of frequency 



P{r) 



p 2 4 
^173 ~~ ~ 



sin(2r + a), p(r) — > for r — > 00 
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contains both monotonous and oscillating constituents thus this equation differs from both 
the oscillator-like equation and the Mathieu-like equation. 

With the time growth the term p 2 /r 4 / 3 drops faster then the absolute value of 
(4/r)sin(2r + a). The condition of their equality gives an estimation of the lower time 
border for the asymptotically long waves: 



jf 



1 



4/3 



T, 



Tic 



n 



ong 



V 



(5.22) 



long " 

From (j5.21j) one can estimate also the upper time border of the asymptotically short waves: 

-.2 



V 



r. 



4/3 



short 



The interval 



1 



At 



Tshort ~ P 



p 6_ p 3/2 



3/2 



(5.23) 



gives an estimation of the transition duration between asymptotic waves. 

We have developed the special asymptotic method which allows us to receive the asymp- 
totic solution of equations like 



y" + 



const 4 

1 sm(2x + a 







with any desired order of accuracy. Formal scheme for construction of general solution is 
similar with one for Mathieu-like equation. At first we don't take into account the oscillating 
correction and solve the oscillator-like equation. Then we add the mixed harmonics into the 
solution, all constants before trigonometrical functions replace by functions in the form of 
power series that depends on n and m, so the whole problem reduces to search of expansion 
coefficients. 

Let us apply this scheme to equation (|5.21j) . We solve the oscillator-like equation 



X 



1 + 



^473 



X 



0. 



Here p 2 /r 4//3 — > steadily in the limit r — > 00, therefore the general asymptotic method is 
applicable for solution of this equation (see Section 2). So we have 



X(r) 



1 



1 + 



27473 



C x e 



3p z 
"SITS 



+ C 2 e 



3p' 
77T73 



Now let's take into account the influence of correction (4/r) sin(2r + a) 
one can write in a form 



2i 



T 



-i(2r+a) 



x = o. 



(5.24) 
Equation (I5~2T1) 
(5.25) 
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According to described scheme the general form of solution in the limit r — > oo is 
1 



X(r) 



27*73 



E 



1 + irirw s=0 L 



, 3p 2 , A ■( 3p 2 . 

A ^ T y(2s+l)(r+a/2) e \^T73 + ^) + B ^ T y~i(2s+l)(r+ a /2) ~ l { ^73 



_ C ( T \ e i(2s+l)(T+a/2) -\-itjz + ^) + jD ^ r ^ e -i(2 S +l)(r+a/2) e i (iff73+^ 



+ 
(5.26) 



where <^i, y?2 = const. In the frequency of equation (|5.25jl there are the terms with powers 
of t multiple to 1/3. Therefore 



fe/3 



k=3 

— oo 



k=3 
— oo 



(5.27) 



k=3 



k=3 



For a function ^^(t) — x/(t — t ) we have obtained the following expression: 
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40m 4m 1/3 
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20m 2 / 3 
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3p 2 
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3p 2 \ ( i 
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4m 2/3 
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(5.28) 
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+ 



i i p 

1 ^ 2m 4 /3(t-t )4/3 



20m 2 / 3 
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t - t 3p 4 (t - t ) 



1/3 



sin 



3p 2 



2m 1 /3(t_t )i/3 



where a is defined in (j2.14j) . Constants C\ and C*2 are found by matching of solution (j5.28j) 
with the short-wave solution ()5.19|) . The matching time is defined from relation ()5.23|) : 
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+ to 



(5.29) 
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Therefore 
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Function of metric perturbations 1)5.28)1 for every regime can be essentially simplified: 
1) Transient regime at t{k) < t < ti ong , where ti ong = p 6 /m + t : 
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2) Long- wave regime at t > ti ong 
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3p 2 



+ C 2 sin 



3p 2 
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Last expression one can also write in the form analogous to (|5.16|) : 
" 3p 2 



40m 
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(5.32) 



. (5.33) 



Obviously the function of long- wave solution (j5.32j) oscillates with the constant amplitude 
as sin (m(t — to) + a /2), however, in contrast to (|5.16J1 . there is a slowly varying modulation 
factor, which approaches to the constant just in the limit t — > oo. 

Theoretical predictions for the parameter of energy density perturbations spectrum are 
of the most interest. This parameter is defined as follows 



A 



01 



5e 
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(5.34) 



The analytical results of calculation of A without the averaging over the initial phases of 
quantum fluctuations are shown in Appendix. 
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VI. CONCLUSION. 



Gauge invariant approach to investigation of linear scalar perturbations of inflaton and 
gravitational fields (SIG-waves) has been developed. We mean the derivation and analytical 
investigation of equation for the single invariant function without resorting any gauge. 
Equation of invariant dynamics (EID) is constructed both in the cosmic and conformal times. 
The transformations of EID coefficients together with background equations for arbitrary 
potential have shown an absence of the mathematical singularities in reality. This fact has 
given a chance for analytical investigations of the EID. This internal property of the invariant 
dynamics is its significant advantage. 

Our approach allows to compare various gauges used by other researchers, and to find 
unambiguous selection criteria of physical and coordinate effects. We have shown that 
the so-called longitudinal gauge commonly used for studying the gravitational instability 
leads to the overestimation of physical effects due to presence of nonphysical proper time 
perturbations in results obtained by using this gauge. 

We have proven that the non-invariant part of energy density perturbations in the syn- 
chronous gauge is equal to zero for any co-ordinates transformations preserving the gauge, 
therefore the calculation of the relative perturbation of energy density 5e/e in the syn- 
chronous gauge gives correct results in contrast to the longitudinal one. 

To a great extent, the technology of EID solving rests upon mathematical symmetry 
properties of invariant dynamics. The general long-wave solution of EID as a functional of 
background solution for arbitrary potential U(<f>) is obtained. We have developed asymp- 
totical methods which allow us to obtain the invariant function of metric perturbations in 
analytical form at all stages of the Universe evolution for arbitrary wave lengths in the 
framework of fixed potential (we use of the simplest model with U(<p) = m 2 (p 2 /2). It is 
important point of our investigation that both the non-perturbed or background charac- 
teristics (Hubble function H(t), inflaton <f>(t) and scalar factor a(t)) and characteristics of 
metric perturbations (invariant functions Jk or ^t) at the inflationary and postinflationary 
stages in the considering model of the Universe are completely defined by parameter 
m and initial values of Hubble function H and scale factor a . Note also that all analytical 
results coincide with numerical solution of EID at various stages and wave numbers with a 
high precision. 

We have also obtained analytical expressions for the energy density perturbations spec- 
trum A(k, t) for all possible wave numbers k and times t. (see Appendix). Amplitude of the 
long-wave spectrum in the case of transition from short waves to long ones occurs at the 
inflationary stage is almost flat, i.e. has the Harrison-Zeldovich form, for arbitrary poten- 
tial U(ip), but there is some tilt that can be compared with one from recent high precision 
data. This result is supposed to be the principal result of our invariant approach. Inflation- 
ary prediction for nearly flat spectrum of density perturbations is in agreement with both 
measurements of the CMB anisotropy and observations of structures in the Universe. 

We are grateful to V.A. Beylin, V.I. Kuksa, O.D. Lalakulich, L.S. Marochnik, A.V. 
Zayakin and O.V. Teryaev for useful discussions and comments. This work is partially 
supported by grant RFBR 03-02-16816. 
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VII. APPENDIX. CALCULATION OF THE ENERGY DENSITY PERTURBA- 
TION SPECTRUM. 



In this section there are the analytical expressions for the parameter of energy density 
perturbations spectrum A(k,t) (|5.34j) : 



A(k,t) 



k 3/2 
V2n 



8 Si 



Expression for relative perturbations of energy density through \l/k has a form: 
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3H 2 



y/2H f y^H 



H 



H 



(7.1) 
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to 



Using the expressions for \l/k at various t and k (|5.14j) . (|5.1fjj) . (|5.19j) . (j5.31j) and (|5.HHj) . we 
have obtained the following expressions for 5ei nv /e and A(k,t): 

1. Short waves at the inflationary stage that is at times < t < t(k), t ^ ti, where 
t(k) is the transition time from short waves to long ones, t± is the time when the 
inflationary stage terminates. In that case we have 
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2. Lon^ waves in the case of the transition from short waves to long ones occurs at the 
inflationary stage that is at times t > t(k), t(k) ^ t\. Here we have 
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(7.6) 



where 



H s = H(t(k)) = J Ho 2 



2m 2 



In 



mao 



is the value of Hubble function at the transition moment. We see that the amplitude 
of this spectrum is almost constant with k for arbitrary potential U(ip). Therefore this 
spectrum has the Harrison- Zeldovich form. Some tilt from flatness contains in H s and 
can be compared with recent high precision data. 
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3. Long waves in the case of the transition from short waves to long ones occurs at the 
postinflationary stage that is at times t > ti ong , t\ ong > t\ : 
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Here instead of k we introduced a new parameter 

p = p ^ = 2 —J^\ » 22/3 ■ 

a mexp J 

Since these perturbations are considered at the postinflationary stage so in the frame- 
work of investigated model 
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Apparently, the long- wave functions 5ei nv /e and A(k,t) are functionals of background 
solution H(t) and depend on time as H(t)/H(t). This statement is correct for any 
potential. 
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4. Short waves at the postinflationary stage that is at times t± < t < t(k), p ^> 2 2 / 3 : 
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Waves of middle lengths (occurs only at the post inflationary stage) that is at times 
h < t(k) < t < ti ong , p > 2 2 / 3 : 
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